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Abstract 

The spinless Falicov-Kimball model is studied by the use of a recently developed projector-based 
renormalization method (PRM) for many-particle Hamiltonians. The method is used to evaluate 
static and dynamic quantities of the one-dimensional model at half-filling. To these belong the 
quasiparticle excitation energy ik and the momentum distribution nj. of the conduction electrons 
and spatial correlation functions of the localized electrons. One of the most remarkable results is 
the appearance of a gap in at the Fermi level of the order of the Coulomb repulsion U, which is 
accompanied by a smooth behavior for n^. The density of states for the conduction electrons and 
the one-particle spectral functions for the localized electrons are also discussed. In both quantities 
a gap opens with increasing U. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.30.Mb 
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I. INTRODUCTION 



The Falicov-Kimball model is a widely used model to study the properties of interacting 
fermionsi and has also been used to describe the valence change transition in YblnCui and 
similar compounds^. (For a recent review see Ref. []). The model considers a lattice of 
localized /-sites, which are either empty or singly-occupied site, and of conduction electron 
sites, which are delocalized by the nearest-neighbor hopping. There are rif /-electrons and 
n c conduction electrons interacting by an on-site Coulomb interaction U. The two types of 
electrons share the common chemical potential fi, which is always adjusted so as to keep 
rif + n c = 1. The local /-charge is a constant of motion but thermal fluctuations can change 
the average /-occupation by transferring electrons or holes from the conduction band to 
/-states. The relative occupation of /-states and conduction states is determined by the 
competition between the entropy of the /- and the conduction states, the excitation energy 
of the f-states, the kinetic energy of band states, and the interaction energy. 

The general solution of the Falicov-Kimball model is not known except in infinite 
dimension^^i. In low dimensions, the zero-temperature phase diagram is highly non-trivial 
and has attracted a lot of theoretical attention (see (^j) for the list of references). In one 
dimension the thermodynamic properties of the model are more or less understood, whereas 
dynamical properties have been much less studied. Kennedy and Lieb^ proved that at low 
enough temperatures the half-filled Falicov-Kimball model for dimensions d > 2 possesses 
long-range order, i.e., the ions form a checkerboard pattern, the same as in the ground state. 
In the ordered phase, the lattice can be divided into two inter-penetrating sublattices A and 
B in such a way that all nearest neighbors of a site from sublattice A belong to sublattice 
B and vice versa. This result holds for any value of the Coulomb interaction U. For one 
dimension longe range order exists only at temperature T = 0. 

In what follows, we first describe the Hamiltonian and the projector-based renormaliza- 
tion method (PRM) which was recently introduced (Sect. I). In Sect. II the method is then 
applied to the one-dimensional Falicov-Kimball model and renormalization equations are de- 
rived for the parameters of the Hamiltonian. The results for static and dynamic quantities 
are discussed in Sect. Ill and IV. 

To simplify the calculation, we consider only the case with translation symmetry. This 
excludes the discussion of physical properties which depend on lattice sites i, i.e. we do 
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not discuss the ordered phase which emerges at low temperatures. However, as discussed in 
Sect. V, we can still evaluate physical quantities which are given by averages over all lattice 
sites. Note that the PRM approach can also be extended to situations with lone ranee order. 
This was shown for instance in Ref. (J8|), where the quantum phase transition for the spinless 
Holstein model was discussed as function of the electron-phonon coupling. 

A. Model 

The spinless Falicov-Kimball model is described by the Hamiltonian 

H = Y U >i c fe + £ f Y Ui + U Y n ^ ' w 

<i,j> i i 

where, c\ and Cj are the creation and annihilation operators for conduction electrons at 
site i, tij are the hopping matrix elements, Ef is the energy level of the localized electrons, 
and rii = f-fi and n\ = cfcj are the occupation number operators of the localized and 
conduction electrons at site i. The summation runs over the N sites of a periodic lattice, 
and U is the interaction strength of the local Coulomb repulsion at sites i between conduction 
and localized electrons. Note that rii at each site i commutes with the Hamiltonian. Thus, 
the /-electron occupation number is a good quantum number, taking only two values 1 or 
0, according to whether or not site i is occupied by the localized /-electron. 

B. Projector-based renormalization method (PRM) 

Let us first discuss a recently introduced many-particle method, the projector-based 
renormalization method (PRM)^, which provides an approximate solution to many-particle 
problems defined by the Hamiltonian H,, 

U = H +Hv (2) 

Here, Hq is the unperturbed part with a known eigenvalue spectrum, Ho\n) = En\n), 
and Tii is the perturbation which does not commute with H.q and has off-diagonal matrix 
elements, (n\Hi\m) ^ 0, but no diagonal elements. In a usual perturbative approach one 
evaluates the corrections to the eigenstates and eigenvalues of TCo due to Hi. The n-th 
order correction requires the evaluation of all matrix elements up to (n|(Hi) n |m), which are 
difficult to calculate. 
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An alternative insight is obtained by making a unitary transformation to a new basis 
which has no transition matrix elements with energy differences larger than some chosen 
cutoff A < A. Here A is the largest energy difference between any two eigenstates of 7Y 
which are connected by H.\. This generates a new Hamiltonian 

H x = e x *He~ x \ (3) 

which can be written as a sum of two terms, 

= 7~Lo,\ + T^i, a, (4) 

such that (m\Hi t \\n) = for \E* — E^\ > A, where E\ and |n) are the new eigenvalues 
and eigenstates of Ho, a- ^i,a is chosen in such a way that it has no diagonal elements with 
respect to 7io,A- The generator of the unitary transformation is denoted by X\ which has 
to be anti-hermitian X\ = —X\. Note that the elimination of high energy transitions may 
also generate new interaction terms which have operator structures different from that of 
Tii. However, they do not connect states with an energy separation larger than A. 

In the PRM method the elimination of transition matrix elements is carried out by 
defining a generalized projection operator Pa which removes from any operator A those 
parts which give rise to 'forbidden' transitions, i.e., 

P\A = ^ \n)(m\{n\A\m), (5) 

m,n 
| | <* 

i.e. we retain in §5§ only the states \m) and |n) with \E* — E^\ < A. The orthogonal 
complement of Pa is Qa = 1 — Pa- Here only dyads \n){m\ contribute with \E n — E m \ > A. 
The generator X\ in ([3]) is determined by the condition 

Qa^a = 0, (6) 

as can be seen by expanding the exponentials in ([3]). This gives 

^{^ + E^ X ^} = °> ( ? ) 

where X A is a superoperator defined by X A 7Y = [X\,TL], X A 2 7Y = [X\, [Xx, H]], etc. In 
lowest order, Q\H\ = reduces to Qa |Wi — LoXl 1 "*! = 0, where L A = [Ho, A] defines 
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the unperturbed Liouville operator L , which commutes with Qa. Obviously, the condition 
Q,\H\ = is satisfied by the expression X^' = Lq 1 Qa7Yi, i.e., to lowest order, xjj^ can 
be obtained from the decomposition of Tii into the eigenmodes of L . The higher order 
correction terms to X^ follow systematically from the higher order commutators of ([7]) (for 
details see, e.g., Ref. |9_|). 

We see from Eq. ([7]) that all transitions in Tii with energy transfers between the original 
cutoff A and the new cutoff A are eliminated in one step. However it is more convenient to 
perform the elimination procedure step-wise, such that each step reduces the cutoff energy 
A by a small amount AA. Thus, the first step removes all the transitions which involve 
energy transfers between (the original cutoff) A and A — AA, the subsequent steps remove 
all transitions larger than A — 2AA, A — 3AA, and so on. The unitary transformation for 
the step from an intermediate cutoff A to the new cutoff A — A A reads (in analogy to ([3])) 

^a-aa = e x ^n x e~ x ^ (8) 

where H\ = H ,x + 'Hi, a and Ha-aa = ^o.a-aa + T^i.A-AA-The generator X Aj aa is now fixed 
by the condition 

Qa-aa^a-aa = 0. (9) 

which specifies that 7^a-aa contains no matrix elements which connect eigenstates of 7Yo,a-aa 
with energy differences larger than A — AA. 

Let us assume that the operator structure of 7i\ is invariant with respect to further 
unitary transformations. Then, we can use Eq. ([7]) to derive difference equations for the 
A-dependence of the coupling constants of the Hamiltonian. These equations connect the 
parameters of the Hamiltonian with cutoff A to those with cutoff A — AA. By using a finite 
number of steps we can proceed to A — > and obtain a set of nonlinear equations for the 
renormalized parameters. The solution determines the final, fully renormalized Hamiltonian 
7Ya^o = Ho.A^cb which depends on the initial parameter values of the original model at 
cutoff A. In this limit all transitions due to the interaction with nonzero transition energies 
have been eliminated so that 7^i,A-»o identically vanishes, i.e. Hi^q = 0. 

The underlying idea of the PRM, namely the elimination of the interaction terms by 



the unitary transformations, has been used before in the literature. For instance, in Ref 
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a unitary transformation has been employed to eliminate off-diagonal matrix elements of 



the Coulomb interaction in order to study ground state properties of the Hubbard model. 
Similarly, in the Schrieffer-Wolff transformation 11 , which maps the Anderson to the Kondo 
model, matrix elements connecting different charge configurations of magnetic ions are elim- 
inated also by use of a unitary transformation. Note however that in the PRM approach 
unitary transformations are performed in small steps in contrast to earlier applications. Note 
also that the approach removes high energy transitions but does not decimate the Hilbert 
space. This is different to the poor man's scaling 12 which removes high energy states. It 
should also be mentioned that the PRM resembles the previous flow equation method™' by 
Wegner and the similarity renormalization approach^ by Glazek and Wilson which can be 
considered as continuous versions of the PRM method. 

II. RENORMALIZATION OF THE FALICOV-KIMBALL MODEL 

In order to apply the PRM we express the Falicov-Kimball interaction in terms of the 
fluctuations with respect to the thermal averages of n\ and n*, and write the diagonal and 
off-diagonal part of H as 

n o = Y1 ( £ f + u i n i)) n i+J2 ( £k + f/ <^» c l Ck - UN ( n i) «)' ( 10 ) 

i k 

and 

Wi = e-« k - k '^ 5 ni 5(clc k ,), (11) 

ikk' 

where bn\ = n\ — (nf), Srii = m — \Tii), (• • •) denotes the thermal averaging with the 
full Hamiltonian H, and c k = (1/y/N) ^e* feBi c] is the Fourier transform of cj. Of course, 
due to the assumed translation invariance (n*) is equivalent to (1/N) $^i( n i)> where iV is 
the number of lattice sites. The sum in (TXT]) is restricted to k ^ k', since the diagonal 
part of the Coulomb repulsion is included in 7i . From translation invariance it follows 
5(c\ck>) = c\ck> — (c k Ck)5k,k'- In principle, the average values in (fTO]) can be defined with 
an arbitrary ensemble since Eq. ([T]) and Eqs. (TTUI) . (TTTT) are equivalent. It has turned out 
that the best choice is to evaluate the averages with respect to the original full Hamiltonian 
H = Ha, which can be done by the procedure explained below. 

The perturbation Hi causes transitions between the eigenstates of Ho by creating 
electron-hole pairs in the conduction band, whereas the number of /-electrons is conserved 
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at each site. Using L Srii 8{c k cy) = (e*. — e^) 5rii 5(c k Ck>), TCi can be decomposed into a 
sum of eigenmodes of the unperturbed Liouville operator. The lowest order solution of the 
generator of the unitary transformation is = Lq 1 Q\Ti,\, which is given by 

tj p i(k—k')Ri 

X X ] = m E - £ «\ - A ) ^ 4 C *'- ( 12 ) 

N lkt> £k ~ £k ' 

Here G(|efc — — A) is the projection operator which removes from Tii all transitions 
with energy transfers larger than A. Note that the same operator form also appears in the 
non-perturbative calculation below. 



A. Renormalized Hamiltonian 

The unitary transformation leads to the renormalized Hamiltonian 7i\ = e Xx 7ie~ Xx , 
which can be written as TC\ = 7io,A + ^i,a, where TCi,\ gives rise to the transitions between 
the eigenstates of Ho,\ but has no diagonal elements. The choice of 7Yo,a is n °t unique, and 
each particular case requires a physical intuition. But in any case it has to be such that the 
thermal averages which appear in the procedure can be evaluated. For the Falicov-Kimball 
model, the choice which emphasizes the weak coupling limit, is 

7<o,\ = £/,aJ] Srn + ^ 9ij,x Sn i Sn j + ^2 6k > x c k Ck + ( 13 ) 

i i^j k 

ikk' 

where Ef t \, gij,\, and Ek,x are the renormalized parameters which are calculated below. Note 
that the interaction parameter U will not change in the renormalization procedure below 
so that a A-dependence of U is not considered from the beginning. The second term in 
( 1131) is a new density- density interaction between the localized electrons which is generated 
during the renormalization procedure. Additional higher order operator terms, which are 
also generated by the renormalization, are neglected. The original model (A = A) is defined 
by the initial condition, according to (fTUl) . ffTTj) 

e f ,K = e f + U{nf) £k,A = e k + U(rii) (15) 

9i h K = E A = -N(m) (e, + U «)). 
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Once we have made the ansatz (fT3j) and (TH1) for the operator structure of 7io,A and 7Yi,a, 
we reduce the cutoff to A — AA to find the new effective Hamiltonian 

Ux-ax = e x ^Hxe~ x ^ = H ,x-ax + Wi,a-aa (16) 

with renormalized parameters. The new Hamiltonian is such that 7^i,a-aa does not give 
rise to transitions with energy transfers larger than (A — AA) (with respect to the new 
unperturbed part Wo,a-aa)- 'Hx-ax should have the same operator structure as H\, however 
with renormalized parameters 

^A-AA = £/,(A-AA) ^ Sni + 9ij,(X-AX) 5Ui Srij + ^ £ft,(A-AA) C ] k C k + -E(A-AA) (17) 
i i^j k 

+ ^^e(A-AA-| £M _ AA -^ A _ AA |)e^- fc '^ (5^5(4^)). 

ikk 1 

B. Unitary transformation 

For the explicit evaluation of the unitary transformation (Tl6!) we need the generator 
X\ ax- We make the following ansatz 

Xx,ax = ±J2 A tt X 5n, 8{c\c k ,). (18) 

k,k',i 

Here 0^f A is a product of 6-functions 

= ©(A-|e M -ev,A|)xe[|e fc>(A -AA)-ejfc/,(A-AA)|-(A-AA)], (19) 

which project all the operators in Xx^ax 'on the shell' between A and A — AA. More precisely, 
the energy differences \£k,x ~ £ k',x\ and |£fc,(A-AA) ~ £ fc',(A-AA) I i n CEH]) re f er to the two different 
Hamiltonians 7i A and Tix-Ax- Therefore, the two B-functions in ([19]) take into account that 
Hx possesses only transition elements with energy transfer |e fcjA — ey ; \\ < A, whereas in 
'Hx-ax no transitions with (^(a-aa) — £ fc',(A-AA)| > A — AA are allowed. The coefficients 
A x '^, x are determined below, using the condition Qa-aa^a-aa = 0. Again, the form of Xx,ax 
is suggested by its lowest order expression according to (|T2l) (compare (|9j)). The unitarity 
of the transformation requires Xx,ax — AA , so that A X '^ X has to be antisymmetric with 
respect to the interchange of k and k! . Expanding the exponentials in (Tl6]) the transformation 
can be written as, 

V 

^a-aa = e x ^Hxe~ x ^ =^ + Y,~\ X a,aa^a (20) 

n=l 
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where X AiA a^a = [Xx t A\,H\], X^ AA 7^ A = [A AiA a, [A Ai aa,^a]], etc. The commutators are 
evaluated in Appendix A. The basic approximation of the PRM approach is to replace some 
of the operators generated by X Ai a A by ensemble averages, so as to keep the structure of H\ 
invariant during the renormalization procedure. This additional factorization enables us to 
resum the series in (1201) . 



C. Renormalization of the coupling constants 

By comparing the respective operators in (|17|) with those obtained from fTSU]) one finds 
renormalization equations for the parameters of T~C\. First, from the prefactors of 5rii one 
obtains 

(l-2(n i )) > p A[AA (1-cos^,) , j x 



+ 



Z7(l — 2(ni)) v^„ A| aa sinw 



A 



V P)A,AA ^ ^ / f \ 



Similarly, by comparing the coefficients of c\ck one finds 

£fe,(A-AA) = £k,\ + ^Y1 Q kf' X ( £ k',x - e k ,x) (1 - cosv^ fc ,) (22) 

The renormalization equation for the interaction parameter <?^ )A between /-electrons is best 
expressed in terms of its Fourier transform g q \, 

, 1 (4 C fc) ~~ (<4+c, C fc+<?) ^A,AA / w-i X \ 

9q,(\-AX) = 9q,X + j 2_> £7~g) l £ fc+9,A ~ ^fc,Aj (1 - COSW fc A . +g J 

+ A= Y {C "° k) '^} Ck+q) - Qii" sin vl k+q (23) 



where (■ • • ) denotes the second and third term on the r.h.s. with the wave vector q replaced 
by q'. In this way, the exact sum rule Yl q 9q,^ = is fulfilled which guarantees that only 
sites with i ^ j contribute to g^. Finally, also the renormalization equation for the energy 
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shift E\ can be obtained which will not be given explicitly. The other new quantities in 
Eqs. f }2~Tj) to fl23|) are defined as follows 



>*» = ^\l qAt ^ 1 A ^ X > CM) = ^E^^H^i)- (24) 



'J 



Finally we have to determine the coefficients A k '^ x of the unitary transformation. They 
follow from the condition Qa-aa^a-aa = 0, which removes from Ha-aa ah the operators 
giving rise to the transitions with energy transfers larger than A — AA. One finds 



n A.AA .A.AA n A,AA 1 N . fn C p (k - k') U \ 

Note that the A kk , x depend on the parameters at cutoff A as well as from the reduced 
cutoff A — AA. The presence of 6^f A on both sides of Eq. fl25|) means that A^ k x has to be 
specified only for those values of k and k! which are 'on the shell' defined by the condition 
G)fcfc^ A = 1. 'Off-the-shell' coefficients can take any value and, in what follows, we use 
A^ = for e^T = 0. Note that expression (I25p includes U to all orders. However, in 
the thermodynamic limit iV — > oo, the coefficients A k ' k , x become linear in U since all higher 
order terms vanish for N —>■ oo. In the actual numerical evaluation of the renormalization 
equations on a lattice of finite size N the excitation energies (£k',x ~ £ k,\) m &y become very 
small for A — ► so that the expansion of A k ' k , x to linear order in U breaks down. In this 



case, the full expression (1251) has to be taken. 

By help of fl25J) one finds for the parameters v kk , of 



sin^, = 2U JC p (k-k>) sign (e M - A ) 



N 



V(?k,x - e k ,, x y + AU 2 C p (k - k')/N 



cosv x kk , = l£k ' X £k '' xl . (26) 

^(e k , x -e k ,, x y + 4U*C p (k-k>)/N 



Equations (1211) to (|2"5j) represent the final renormalization equations for the parameters of 
the Hamiltonian as the energy cutoff is reduced from A to A— AA. The overall renormalization 
starts from the original cutoff A = A, where the initial parameters are given by ( 1T51) . For the 
parameters smv kk ,, cosv kk ,, and C p {k — k'\ which contain expectation values, an appropriate 
choice is taken at cutoff A. Using these quantities we obtain from the renormalization 
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equations (12T]) to (|23|) the renormalized parameters at the smaller cutoff A — AA. The 
procedure is repeated until the cutoff is reduced to zero. At A = 0, we obtain the completely 
renormalized Hamiltonian H : = 7^(a=o), which describes an effectively decoupled system of 
free conduction electrons and interacting localized electrons, 

H = Hc + Hf + E, (27) 

where 

H c = ^£fc c\c k , and H f = if^Srii + ^gij SriiSrij. (28) 

k i ij^j 

The fully renormalized parameters are defined as = £k,(\=o),£f — e f,(\=o)i9ij — 9ij,{\=o) 
and E = E^\=o)- Note that for A —>■ the interaction 7ii 5 A has completely vanished due 
to the ©-function in ( TH1) . Using H we next can recalculate the statistical averages (cj,Cfc), 
C p (k — k'), etc., which appear in Eqs. (T2T1) - (1231) . to find new values for the renormalized 
coupling parameters. This procedure is repeated until the self-consistency is reached. Note 
that the Hamiltonian Tif for the localized electrons is usually denoted as the lattice gas 
model. Due to the dependence of on i — j, the interaction in Hf is not necessarily 
restricted to nearest neighbors. 

Let us add two side remarks: The first one concerns the validity of the present approach. 
As is obvious from the renormalization procedure in Sect. II we are here dealing with a 
weak-coupling approach. The Coulomb interaction U is taken as the 'small' quantity which 
is successively integrated out in the procedure. However in each renormalization step not 
only the lowest order in U is included. Instead, certain renormalization contributions up 
to infinite order in U are taken into account. Therefore, we expect the present approach 
to be valid for U values which are less or at most of the order of the hopping term t. For 
larger values of U either newly generated renormalization terms would have to be added. 
Or, more appropriate for large U would be to take the dominant U term as part of a new 
'unperturbed' Hamiltonian 7io and the hopping term as part of the new perturbation Tii. 

The second remark applies to the optimal choice of the renormalization interval AA. It 
is best determined by the size of the system, i.e., by the requirement that the 'energy shell' 
contains only a few fc-points. Thus, the number of AA intervals should be of the order of 
the lattice size N. For a larger number of intervals the result would be the same since more 
steps occur without renormalization since no transition energies fit into the interval. On the 
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other hand a smaller number of AA intervals leads to a quite large renormalization in some 
steps since many fc-points could contribute in these intervals. In practise, for a given iV we 
have increased the number of intervals until the renormalization does not change any more. 



III. EVALUATION OF STATIC EXPECTATION VALUES 

The renormalization equations depend on expectation values (rij), (niiij) and (c k Ck) 
which have been defined with the full Hamiltonian TC. Since thermal averages are invariant 
with respect to a unitary transformations, we can do the averaging with the renormalized 
Hamiltonian 7Y, provided we also transform the operators, 

(A) = (A(X)) m = (A(X - 0))fl, A(X) = e x *Ae~ x \ (29) 



To evaluate (c\,Ck) one needs the renormalization equations for the operators (c|.c&)(A) = 
c\(X) Cfc(A), where Cfc(A) = e Xx - AX Cke~ Xx < Ax . For Cfc(A) we use the following ansatz, 

Cfc(A) = a KX c k + — ^Pk,k>,\ e i{k - k ' )Rl 5n iCk ,, (30) 

k',i 

where a^x and Pk,k',X are A-dependent coefficients. The operator structure of c&(A) in f[3"Uj) 
is again taken from the first order expansion in U. All higher order terms are neglected. 
The renormalization equations for the new parameters can be derived in close analogy to 
the renormalization of Tlx- One finds, 



a k ,(x-AX) = «fc,A - a k ,\ cos (A x k £ x \j Cp ^ N ^ )] 6^t A (31) 



fe' 



and 



/ ^ • / aX.AX C p (k — fc'h-, ~A,AA 
Pkk',(X—AX) = Pkk',X ~ ttfc.A \/ (n> (fc _ fc/) Sm (A** V ]V ^ fcfc ' 



- (1 - cos [A^WyQ] ) Qt>\ (32) 

where A x 'f x and C p {k — k') are defined by Eq. (|25|) and (124)) . respectively. The initial values 
are again those of the original model (A = A), a k) \ = 1, P kk ',A = 0. Eqs. fl30|) - fl32l determine 
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the transformation behavior of Cfc(A). By the use of (!30|) we finally obtain 

n k = (c[c k ) = |a fc | 2 (4.c fc )^ + — \Pk,k'\ 2c p( k - k ')( c l' c k')n c ^ ( 33 ) 

k> 

where a k = otfc ) (A=o) an d /3fc,fe' = Pkk',(\=o), an d ( c I c fe)^ is the Fermi function formed with 
We, <4c fc >^ = (e^ + I)" 1 . 

The occupation number operators rii of the localized electrons commute with X\. Thus 
rii will not be renormalized, rij(A) = rij, (rij?ij)(A) = n^n.,-, i.e., 

% = ( n i) = ( n i)n f (driiSrij) = (dniSnj)^ . (34) 

Here the expectation values are formed with the lattice gas model Tif since the conduction 
and the localized part of the renormalized Hamiltonian 7i are decoupled. The expectation 
values in ([541) can be found by direct numerical evaluation in the local occupation number 
representation. 

IV. DYNAMIC CORRELATION FUNCTIONS 

Dynamical quantities can be evaluated by the PRM approach along the same lines. In 
what follows we shall discuss the spectral functions of the conduction and localized electrons. 
The one-particle spectral functions of conduction electrons A^ (oo) and (u) are defined as 

At{u) = — / (c k (t) 4) e**dt , A-{u) = — / (4 c k {t)) e**dt , (35) 

•J — oo */ — oo 

where A\{io) describes the creation of an electron k at time zero and its annihilation at time 
t whereas in A^(u) first an electron is annihilated. It is well-known that A^(u) and A~j~(u) 
can be measured by inverse photoemission (IPE) and by photoemission (PE) experiments. 
To evaluate A^(u) and A k (u) we again exploit the property (129|) of expectation values and 
the ansatz fl30l) for c/.(A) and 4(A). We find 

^* H = |« fc | 2 (1 - (4 c *)« e ) % - g fc) + Tr Yl \M* - - ( c ^)^) % - e*,) 

fc' 

(36) 

A£ (w) = |a fc | 2 <4 c fe)^ c - £fc) + ^ |/3 fc fc'| 2 C P (A; - k') (4' c A;')w c <K W ~ 

fc' 
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where (c' k ,Ck')^ is the Fermi distribution. Note that in fl36|) we have already reduced 
the cutoff to A — ► 0. The first term in A^(u) and A^{uj) describes the coherent one- 
electron excitation which is the only excitation in a free electron gas model. The remaining 
contributions are incoherent excitations due to the coupling of c k to different conduction 
electrons with wave vectors k'. The coupling is mediated by density fluctuations of local 
electrons which lead to the appearance of the charge susceptibility C p (k — k') of the localized 
electrons in (1361) . 

The spectral sum rule can be found by integrating (A^(u) + A k {oj)) over u and using 
[4> c *]+ = !> 

du (A+(u) + AM) = \a k \ 2 + - \Pkk>\ 2 C p (k - k'). (37) 
■°° k' 

This relation can also be derived directly from (1301) . 

Next, we evaluate the averaged one-particle spectral functions B + {uj) and B~(u) of lo- 
calized electrons, 

B» = — J (/<(*) f}) e**dt B-(u) = —J Jf} /,(*)) e **dt. (38) 

Since the local /-charge is a constant of motion, all correlation functions between different 
sites i and j vanish. For the evaluation we again use (I29I) and the following ansatz for the 
transformed /-operator 

/i(A) = vxh + ^ $> W e^-V* Sfa^fi, (39) 

kk' 

where v\ and [i k ,k',\ are the A-dependent coefficients. As before, the operator structure of 
the second term is assumed to agree with the one obtained from the first order contribution 
in U and higher order terms have been neglected. The initial parameter values are (A = A): 
z/a = IjZ-ifcfc'.A = 0. The renormalization equations for v\ and fik,k',x can easily be derived. 
In the limit A — > we find the following result 

B + (u) = v 2 B+{u) + l^w| 2 /(4)(l - f(e k> ))B+(uj + e k - e k ') (40) 

kk' 

B~(u) = u 2 B^{u) + Yl lml 2 /(e*)(l - f(h'))B (u - (e k - £»). (41) 

kk' 
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We have introduced new correlation functions 

1 poo -| poo 

£o» = ^ J Jte™ (fi(t)fl) nj B (u) = — J Jte™ (fUi(t))n f . 

(42) 

which are defined with the lattice gas model TCf. In Appendix B the functions Bq(lo) are 
evaluated by use of the retarded one-/ electron Green function Qq{oj), defined with 7if, 

/oo 
dte i ^ t ([f i (t),f!} + )n f (V = + ). (43) 

V. RESULTS 

A. Static properties 

We now discuss the results obtained by numerical evaluation of the renormalization equa- 
tions for the one- dimensional lattice. Note that in the following all energy values are given 
in units of the hopping matrix element t. One first has to calculate the fully renormalized 
parameters £f, and g^ which determine the effective Hamiltonian 7i = 7i c + TCf + E of 
( I28p . All static expectation values which enter the renormalization equations are evaluated 
self-consistently with the full Hamiltonian, following the procedure outlined in Sections III 
and IV. The effective Hamiltonian Tif describes the lattice gas model of ions which interact 
via a possibly long range interaction g^. Thus, we use an additional Monte Carlo simulation 
to evaluate the expectation values rJ = (rii)^ and C//(«) = (SriiSrij) = (SriiSrij)^, which 
enter the renormalization equations together with \c\ck) _ = [1 + e^ £fc ] _1 . 

As already mentioned, we start by making a guess for the initial values of n/, C p (k — k') 
and (c\.Ck), and calculate the renormalized parameters for A — > 0. That is, we find 
Sf = £/,a=o> £k = £k,\=o, and g q = g q> x=o by solving numerically the renormalization equa- 
tions (12~TT)-(|23"1). Next, using the effective coupling g q we find the new values of the spatial 
interactions #y and then recalculate the localized electron expectation values rif and C//(k) 
using a classical Monte Carlo simulation. The Fourier transform of Cff(n) gives an approx- 
imate expression for C p (k — k'). Using new values of rif, C p (k — k') and (c\.Ck) we iterate 
( |2T|) - (|23|) until a self-consistent solution is reached. 

In the following we consider only the one-dimensional Falicov-Kimball model at half- 
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filling, 

n/ + n c =l, (44) 

where rif = (n*) and n c = ^/N^2 k (c\ck) is the averaged local occupation numbers of 
conduction electrons. Condition (1441) represents a point of special interest for valence and 
possibly metal-insulator transitions, caused by the promotion of electrons from localized / 
orbitals ( f n — > f n+1 ) to the conduction band states 15 . 

1. Temperature dependence 

In Fig. [I] the averaged local /-occupation number rif and the renormalized /-energy if 
from the PRM approach are shown as function of temperature for U = 1.0 and several values 
of unrenormalized /-level position Sf. Note, the chemical potential \x has to be adjusted at 
each temperature T, so as to fulfill condition ( 1441) . In the evaluation procedure this requires 
extra numerical cycles, in addition to those from the self-consistent determination of the 
expectation values. Therefore, only a rather small system of 20 lattice sites has been used 
to generate Fig. [TJ Also the considerably increased effort in the MC sampling for the lattice 
gas part in 7io,A a t higher temperatures suggest to restrict oneself to a rather small system. 
This will be discussed below in some more detail. As expected, for small values of the bare 
/-level Sf (dashed curve in Fig. CEK), the average /-occupation rif increases with decreasing 
temperature and tends to 1 at T = 0. The corresponding values of the renormalized /-level 
remain below the chemical potential (if < in Fig. [lb), so that all the sites are occupied 
with one /-electron. In contrast, for large values of Sf (solid curves) rif tends to zero for 
T — > and if stays above the chemical potential (if > 0). Thus all the /-sites are empty 
and the renormalized conduction electron band is completely filled. Finally, at intermediate 
values of the bare £/- value (dashed dotted curves) the / occupation at T = tends to some 
value between and 1, signalizing an intermediate /-valence. In that case, the renormalized 
/-level approaches the chemical potential (if = 0) at T = 0. Below the chemical potential all 
the conduction electron states are occupied and /-states only partially, due to the condition 
rif + n c = 1. 
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FIG. 1: (a) /-electron occupation number rif as function of temperature for U = 1.0 and for three 
different values of the bare /-level £/ = 0.5 (dashed), ej = 1.2 (dashed dotted), and £/ = 3.5 
(solid). The open circles (sf = 0.5) and open squares {et = 1.2) refer to ED results of reference^. 
Note the good agreement with the PRM results. The values of £f are measured from the lower 
edge of the conduction band, (b) Renormalized /-electron energy ej as function of temperature 
for the same values of U and £y-values as in panel (a). 

2. Low temperature properties at fillings n/ = 1/2 and rif = 1/3 

Next, we discuss static properties for rif = 1/2 and 1/3 at low temperatures. The 
parameters \x and Ej have to be adjusted in such a way that the condition (I44p is fulfilled for 
a given value of rif. Note that in a Monte Carlo simulation a nonzero temperature has to be 
used, which is set equal to k&T = 0.1. Also the most dominant ground-state configurations 
have been incorporated in the Monte-Carlo simulation so that the system size could be 
extended to N = 150 sites. For rif = 1/2 the chemical potential \i is in the middle of the 
renormalized conduction band 

In Fig. [2b the result for obtained from the PRM is plotted as function of k for n c — 1/2 
and for different values of the Coulomb repulsion U. The results for n& from the PRM 
approach are very similar to those from Ref.— (for the same values of U). 

In Fig. [2^, the renormalized one-particle energy e^ is shown as function of momentum k 
for n c = 1/2 (rif = 1/2), and for the same interaction values U as in Fig. [2b. Note, that 
can be considered as the excitation energy of a quasiparticle of the original interacting model 
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FIG. 2: (a) Renormalized conduction electron energy for different Coulomb repulsions U (= 
0.1, 1.0,2.0,3.0) and n c = rif = 1/2. Note the appearance of a gap of magnitude U at the Fermi 
momentum kp = tt/2. (b) Same quantity for n c = 2/3 (rif = 1/3) for U = 0.5, 1.0,2.0,2.5. 
shows two gaps where the one at the Fermi momentum kp = (2/3)-7r is now smaller than U .(c) 
Momentum distribution = (ctc^) for the same parameters as in panel (a) and rif = n c = 1/2. 
The smooth behavior of n& at the Fermi momentum follows from the gap in (d) for two 
different temperatures T = 0.1 (solid line) and T = 0.2 (dashed dotted line) and U = 2.0, nf = 1/2. 
Note that the gap at kp for T = 0.1 is completely smeared out for the higher temperature T = 0.2. 

([1]). This follows from the diagonal form of 7i c in (|28|) and the fact that 7i is obtained from 
the original model by a unitary transformation. The main feature of Sk is the appearance 
of a gap at the Fermi level, which opens for any, even small, value of the repulsion energy 
U. Such a gap was recently anticipated by Farkasovsky^ from the correspondence of the 
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Falicov-Kimball model to the Tomonaga-Luttinger fermions. There, the smooth behavior 
of rik with k also gives rise to a gap in the charge excitations. Note that for n c = 2/3 
(rif = 1/3) (Fig. [2b) again a gap opens at the Fermi energy for hp = (2/3)7r for all values of 
U. Thus, the Luttinger theorem is fulfilled even for the present case of an insulating state. 
By closer inspection of Figs. [2h.,b one finds that the gaps at kp scale linearly with U both 
for n c = 1/2 and n c = 2/3. The proportionality constant is almost 1 for n c = 1/2 whereas 
for n c = 2/3 it is somewhat smaller than 1. 




FIG. 3: (a) Correlation function Cff{n) = (<5nj<5nj+ K ) of localized electrons for rif = 1/2. For 
rif = 1/3 we obtain an analogous result. The length k characterizes the distance between the 
sites. The results correspond to averages (1/iV) ]T^(<5nj5ni+ K ) over all lattice sites (see text), (b) 
Correlation function Cjj(k) = (5rii5ni +K ) of localized electrons for a chain with 20 lattice sites for 
three different temperatures T = 0.20, 0.13, 0.07 and U = 1.0 and n f = 1/2. 

Fig. [3^ shows the spatial correlation functions Cff(n) = {5n,i5n i+K ) between the localized 
electrons for the same parameters as in Fig. [2h. The index k denotes the distance between 
the two involved lattice sites. Note that for nf = 1/2 the values of the correlation function 
are always the same at every second site (k = 2, 4, 6, • ■ ■ ) whereas for all sites in between 
(k = 1, 3, 5, • • • ) they differ only by a minus sign. 
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3. Discussion 

Note that the long-range correlations develop already for arbitrarily small values of U. 
This behavior should be related to the long range order of the one-dimensional Falicov- 
Kimball model at zero temperature 6 . Note that in the present calculations the temperature 
is chosen small but finite {ksT = 0.1), so that in the strict sense no infinitely long range 
order exists. The origin of the long-range correlations can easily be understood from a low 
order U expansion of the effective /-coupling. For the Fourier transform g q of the local 
quantity g^ one finds by use of ff23l) . ff24|) and (125]) 



where the last term excludes equal sites i = j in gij. As follows from (l4"5l) . the interaction 
between the localized electrons is caused from the coupling to particle-hole excitations 
e k+q,\ ~ £fc,A where the wavevector runs over the whole Brillouin zone. Here, the energies Ek,\ 
and £k+q,\ have to be either below or above the Fermi level, i.e. \k\ < kp and \k + q\ > kp 
or \k\ > kp and \k + q\ < kp. Note that for q = ir and rif = 1/2 this condition is fulfilled for 
all /c-values of the Brillouin zone so that the most dominant renormalization of g q occurs for 
q = tt. In contrast, for g-values different from 7r always less k terms contribute. For instance, 
for q ph only k points from the sum in (T4"5l) can contribute which are located in a small 
region around the Fermi momentum kp. As a result of the dominant coupling g q at q = it 
one is immediately lead to a infinitely ranged spatial correlation between the /-electrons 
having an alternating sign from site to site. This explains the long-range behavior in 
Fig. Note again that this feature is not based on a strong-coupling argument. Instead, 
the long-ranged correlations are observed already for small values for U . 

The alternating behavior of Cff(n) = (5ni5n i+K } with n can be interpreted as follows: 
At zero temperatures the model has two equivalent degenerate ground states of period two 
(for rif = 1/2), separated by a potential barrier, with the arrays {1010. ..10} and {0101. ..01} 
of the localized ions. Therefore, we might use the following ansatz for the state vector of 




(45) 
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the system at very low but finite temperatures, 

|$( 2 )) = -L(|10101...10) + |0101...01)) . (46) 
v2 

In order to obey translation invariance the sum of the two generated ground states have 
been taken and ground-state fluctuations have been neglected. Note that the two states are 
never connected by nonvanishing matrix elements when physical quantities are evaluated. To 
proceed, we use (IIS]) to evaluate the correlation function (nifii +K ) of Fig. [31 As explained, the 
result corresponds to the average (1/N) J2i( n i n i+^) over a U lattice sites and is independent 
of i. From ansatz (j4"6l) one finds 

f i : k= 2,4,6,... , 

(uin i+K ) = { 2 (47) 

[ : k = 1,3,5,... 
With (rii) = 1/2 one obtains for C//(«) = (Srii 5n i+K ) 

[ \ : k = 2,4,6,... 
Cff(K) = I I _ (48) 
I — 4 : ^ — 1, 3, 5, ... 

Obviously the result f|48|) agrees very well with the outcome from the PRM approach ac- 
cording to Fig. [3^ for small but also for rather large values of U. A similar reasoning may 
also be given for n/ = 1/3 (n c = 2/3). Here, the degenerate ground-states have period three 
with arrays {100100.. .100}, {010010.. .010}, and {001001. ..001} in analogy to (136). 

Note that the amplitude of the /-correlation function (5rii 5n i+K ) is almost independent 
of U. Already for small [/-values the /-electrons are ordered according to the qualitative 
description of Fig. [3k. The behavior of the conduction electrons is quite different. The cor- 
relation function C c f(n) = (^n^n i+K ^ describes the spatial correlation between a conduction 
electron at site i and an / electron at site i + k. An evaluation for C c /(«;) correct up to first 
order in U gives 

C ^ = ^Ei^i 2 < c ^>^ + ^(- 1 ) K E«^< c ^k+°( f/2 ) ( 49 ) 

k k 

(Q — 7r), where the renormalized parameter /3k+Q ~ U is negative (compare (|32[) ). Obvi- 
ously, the second term in (j4*9l) describes an oscillating behavior with k of c/-correlations 
which favour the presence of conduction electrons in between the / sites. However, in 
contrast to the oscillating behavior of C//(/«), the amplitude of this charge-density wave-like 
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behavior of the correlation function is proportional to U. 

For an interpretation of the temperature dependence of Cff(n) note that the two arrays 
{1010. ..10} and {0101. ..01} (for rif = 1/2) of the localized ions are equivalent to the two 
Neel states | Tit! " ' ' t!) an d | itit ' ' ' it) °f the one-dimensional Ising antiferromagnet. 
For the latter model the spatial spin correlation function (S?S? +K ) decays exponentially at 
low temperatures as exp (— k/£) due to excitation of kinks or domain walls. The correlation 
length is given by £ ~ exp (/3\ J z \/2), which shows that it becomes infinite only exactly at 
zero temperature. J z /2 is the characteristic energy for the creation of domain walls^. For 
the one-dimensional Falicov-Kimball model a quite similar behavior is expected for the 
correlation function Cff(n). Here, also domain wall states can be thermally excited. They 
are obtained by interchanging the positions of neighboring /- and conduction electrons. 
In order to investigate the temperature dependence of the correlation length we have 
evaluated C//(k) = (5ni5n i+K ) for different temperatures T (Fig. [3b). Note that for higher 
temperatures excited states would have to be included in the MC-sampling which leads to 
considerable numerical increase. Therefore, instead of taking a system with 150 lattice sites 
as before a rather small system with only 20 sites was investigated where an exact diago- 
nalization was performed. Fig. [3b shows that for the two higher temperatures T = 0.20 and 
T = 0.13 the correlation function decays roughly exponentially to zero within a few lattice 
sites whereas for the lowest temperature T = 0.07 no decay is observed. Thus, in the latter 
case the correlation length £ obviously exceeds the length of the lattice. The same feature 
was already found in Fig. [3h, where no decay was observed up to K ma x ~ 70. We can con- 
clude that for the lowest temperature no excited states contribute to the correlation function 
for the given lattice sizes. Note that this result is not a consequence of approximations used 
in the present PRM approach. Instead at the lowest temperature the correlation length 
is larger than the lattice extension but should be still finite. Only exactly at T = the 
correlation length £ should be infinite. An equivalent behavior is also found for the itinerant 
electron-dispersion ik as function of temperature. In Fig. [2H the dispersion of ik is shown 
for two temperatures T = 0.1 and T = 0.2 for U = 2.0 and rif = 1/2. Whereas at T = 0.1 
the gap at kp is clearly seen it is interesting to note that it is already smeared out for T = 0.2. 

Note that an infinitely large correlation length £ leads for the Fourier transform C p (q) 
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of Cff(n) to C p (q) ~ N5 qj7r . Instead, for a finite correlation length £ one would find a 
'softened' Kronneker function around q = ir within a range of the order of the inverse 
correlation length Thus the gap in the one-particle excitation energy of Figs. [2^L,b 
also softens over a range in /c-space of the same order It folllows that the system is 
metallic at arbitrarily small temperature. The softening of the gap could hardly be seen 
in Figs. [2k,b for a correlation length larger than the lattice extension of, may be, 100 
times the lattice constant. Thus, the results for Ek and other quantities are extremely good 
approximations also for the T = limit. 



B. Dynamical properties at low T 

Next, we evaluate the dynamical properties using the PRM approach. First we consider 
the electronic density of states (DOS) of conduction electrons 

1 1 1 f°° 

Pc(w) = jy£(4t(")+^(")) = ]v£^ / ([c k (t),cl] + )e^dt. (50) 
k k 

where A\ (u) and (u) were defined in fl35l) . The results from the PRM approach for 
the neutral case rif = n c = 1/2 are shown in Figs. and Hfc for three different values 
of U. For small U the DOS resembles that of free conduction electrons in one dimension 
with a divergence at the lower and upper band edge (dotted line of Fig. [Da,). For larger U 
the opening of a charge gap of order U around the Fermi level at uj = can be observed 
(Fig. and Fig. 2b). For still larger U only charge excitations away from the Fermi level 
are left. This behavior can easily be understood from the former expressions ( |36|) for the two 
spectral functions. Note that for rather small values of U the first terms in both equations 
are the dominant ones since the prefactors |/3kk'| 2 of the second terms are of order U 2 . 
From the gap in the renormalized quasiparticle energy Ek of Fig. [2^l the occurrence of a 
gap in p c (w) immediately follows. Moreover, since the dispersion of around the Fermi 
momentum flattens with increasing U also the peaks in the DOS at the edges of the gap can 
be explained. The result for the larger U value in Fig. 0b show a similar behavior as for the 
lower U's. Note however, for these cases our weak-coupling approach might be at the limit 
of its validity. 

In Fig. Hh also the exact result for the DOS from Ref.— at temperature T = is shown 
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FIG. 4: (a) DOS p c {u) as function of u for U = 0.1 (dotted line) and U = 1.0 (solid line) and 
n>f = n c = 1/2. Note the opening of a gap in the excitation spectrum at the Fermi level u = 0. The 
dashed line is the exact result for T = as taken from Ref.—. (b) DOS p c (oj) for a different filling 
nf = 1/3 (n c = 2/3) for two values J7 = 1.0 (solid line) and U = 2.0 (dashed line). There are now 
two gaps found in the charge spectrum, (c) Same quantity as in panel (a) with nj = n c = 1/2 for 
a larger U value U = 2.0 (solid line). Besides the low temperature result the DOS is also shown for 
a higher temperature T = 0.2 (dashed dotted line). Note that the low-temperature gap at uj = 
is smeared out for T = 0.2. 

(dashed line). There, the divergencies at the band edges appear since an infinitely large 
system size was considered in that approach. In contrast the renormalization equations in 
the PRM approach were solved on a finite lattice by which the divergencies become smooth. 
On the other hand, in order to satisfy the sum rule the spectral function within the bands is 
enlarged in the PRM result. In a recent Monte- Carlo study^ the two-dimensional Falicov- 
Kimball model for the same concentrations n/ = n c = 1/2 was discussed. It shows quite 
similar features for the DOS as our PRM result for dimension d = 1. In addition in Ref.— 
also temperature-dependent effects have been discussed. 

In Fig. HJd the DOS p c (u}) is shown for rif = 1/3 (n c = 2/3) for our lowest temperature 
T = 0.1. As expected from the previous results for in Fig. [2Jd, there are now two gaps 
found in the charge excitation spectrum, one of which is again located at the Fermi level. 
Note that in Fig. H)d the density of states is symmetric with respect to the band center. This 
symmetry can easily be understood as follows: First, at low temperatures the correlation 
function (5ni5rii +t ^ is invariant with respect to translations by three lattice sites (compare 
the lower panel of Fig. [3]). Therefore, C p {k — k') ~ iV 8k-k',±2-K/z at l° w enough temperatures. 
By exploiting this relation and Eq. (j22j) one finds that all k points which are symmetric to 
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±7r/2 lead to the same renormalization of the quasiparticle energy e^x but with an opposite 
sign. Thus, the renormalized quasiparticle energy e& becomes antisymmetric, since also the 
original was antisymmetric relative to ±7r/2 and the DOS becomes symmetric relative 
to the band center in the present approach. This result is in contrast to the expectation 
that for occupation n c away from half-filling the DOS should be asymmetric. In fact, one 
can assure oneself that higher order terms to Xa,aa would give rise to a small additional 
asymmetric contribution to the DOS in our approach. Note that such a small asymmetry is 
found in an exact result 2 ^- for the DOS in one dimension at T = for n c = 2/3. 

Finally, in Fig. 0b also the influence of temperature on p c {oj) is shown for n/ = n c — 1/2. 
Note that at T = 0.1 the gap in the DOS is clearly seen. However, by increasing the 
temperature by a factor of 2 the gap becomes completely smeared out. 
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FIG. 5: One-particle spectral functions A£ (u>) (solid line) and A k (u>) (dashed line) at low temper- 
ature (T = 0.1) as function of u for different k vectors, (a) k > kp, (b) k > kp, and (c) k < kp. 
(U = 1.0, n f = n c = 1/2.) 

In Figs. Eh-c the one-particle spectral functions A^(u) and Aj~(u) for different k values 
are shown for the lowest possible temperature of T = 0.1 and Uf — n c — 1/2. For k > kp 
(Fig. [5^) a coherent excitation is found in A^(u) and a small incoherent excitation in A k {uo). 
To understand the latter feature, note that only a single k' value, k! = k — tc, contributes 
to the incoherent part in (1361) due to C p (k — k') ~ NS^-w^- The same remains true when k 
comes closer to kp (Fig. [Sb). In this case, however, the incoherent contribution from A^(u) 
has gained weight. The reason is the renormalization factor | | 2 in equation (1361) . It 
becomes stronger when k and also k' approaches the Fermi momentum as follows from (|32|) 
and (j2"5j) . Finally, for k < kp (Fig. [5b) a coherent contribution is found in A k {u) and an 
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incoherent excitation in A^(u). 

In Figs. [6^-c the one-particle spectral functions B + (u) (black) and B~(oj) (red) of the 
localized electrons averaged over all lattice sites from section IV are shown for different 
values of U and Uf — n c — 1/2. For the two lower [/-values (U = 0.1, 1.0) at a rather low 
temperature T = 0.1 the coherent parts (of weight z> 2 ) dominate in Fig. Eh the spectrum 
(compare (139]) ). In contrast, for the two larger [/-values (U = 2.0,3.0) the incoherent 
excitations have gained considerable weight (Fig. [6]d and Fig. [6b). Note that at a higher 
temperature T = 0.2 the incoherent part is smeared out over a wide frequency range for U = 
2.0 (dashed dotted curve in Fig. [6b). The same feature would also be expected for U = 3.0. 
However, no stable solution of the renormalization equations was found for this temperature. 
The discussed behavior can again be understood from the former PRM expressions (|40p . (|4ip 
for the spectral functions B + (lv) and B~(u). The first term in both equations describes the 
coherent excitation of weight v 2 at u ~ if = whereas the second term follows from the 
coupling of fi to electronic particle- hole excitations (compare ( 139]) ). Their weight jl\ k , is at 
least of order U 2 . The gaps in Figs. Eh-c of order U immediately follow from the gap in the 
quasiparticle energy 

The same quantities B + (uj) (black) and B~(uj) (red) are also shown in Fig. [7] for /-filling 
Uf = 1/3 (n c = 2/3) and U = 1.0,2.0. The gaps between the respective coherent and 
incoherent parts again correspond to the gap of the quasiparticle excitation in Fig. [2)3 at the 
Fermi level. Note that the spectrum is no longer symmetric around uj = 0. The different 
weights of B + (u) and B~(u) in Fig. [7] follow from the Fermi functions in ( 140]) and ( f4~T]) . 

We have also checked the /-sum rule 



for the spectral functions in Figs. [6] and [7J We found that for the two small U values in 
Fig. [6^ the sum rule was perfectly fulfilled whereas for the spectral functions with U = 2.0 
and 3.0 (Fig. [6Jd and Fig. [6fc) the agreement was less satisfactory. For instance, for U = 3.0 
a value of 0.9 instead of 1 was found. The origin of these deviations can be understood as 
follows. First, from Eqs. ( 140]) and ( f4~T]) one finds that ( |5T]) is fulfilled if 



kk' 

is satisfied for all values of A. For the original model relation (152]) is trivially valid. By 
restricting oneself to the lowest order processes in the renormalization equations for the 




(51) 




(52) 
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A-dependent coefficients in ( 139!) one can show that the sum rule is fulfilled for (A — AA) 
if it was fulfilled for A. However, there are deviations due to higher order renormalization 
contributions for larger U and small values of A as follows from expression (I25p for A^' k , x . 
Also important for the deviations is the chosen simple operator structure of the generator 
Xx t A\ from fll2p . Note however that higher order processes in U are suppressed for larger 
system sizes N. Thus, in order to control the influence of higher order terms a finite size 
scaling in N should be done. Instead, in the present approach a simpler truncation procedure 
has been used. The renormalization was stopped at the smallest A value at which the 
influence of higher order terms was still negligible. This leads to the observed difference in 
the sum rule as compared to the exact result. 
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FIG. 6: (a) One-particle spectral functions B + (uj) (black) and B~(lo) (red) of /-electrons as 
function of u for U = 0.1 (dashed line) and U = 1.0 (solid line), and nj = n c = 1/2, T = 0.1. 
For both U values the coherent excitation dominates the spectrum. (b),(c) Same quantity as in 
Fig. [6ti for two larger values of U, U = 2.0 (b) and U = 3.0 (c) and T = 0.1 (solid lines). Note that 
now the incoherent excitations have gained considerable weight as compared to the case of panel 
(a). For U = 2.0 (b) the spectral function is also shown for a higher temperature T = 0.2 (dashed 
dotted line). 



VI. CONCLUSIONS 

In summary, in this paper we have applied a recently developed projector-based renor- 
malization method (PRM) to the spinless Falicov-Kimball model in one dimension and at 
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FIG. 7: Same quantity as in Fig. [6]for U = 1.0 (a) and U = 2.0 (b) and nj = 1/3. 



half-filling nf + n c — 1. As a starting point in the formalism, we have taken the one-particle 
part of the initial Hamiltonian (TjQ) as unperturbed Hamiltonian Tio and the local Coulomb 
interaction U between conduction and localized electrons as interaction 7i\. The elimination 
of the latter part leads to the renormalization of the Hamiltonian. Thereby a new density- 
density interaction between the localized electrons is generated. Finally, a number of static 
and dynamic quantities has been calculated. Note that the present formuation of the PRM 
is a weak-coupling approach which however should lead to reliable results for ?7-values up to 
the order of t. It is shown that the investigation of the small temperature behavior enables 
us also to extract valuable information about the zero temperature properties of static and 
dynamic quantities. 
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APPENDIX A: OPERATOR TRANSFORMATION OF H\ 



In this appendix we evaluate the renormalized Falicov-Kimball Hamiltonian obtained 
by the unitary transformation W(a-aa) = ex P{-^A,AA} 7~lx exp{— X^aaK where X Ai aa is 
given by (1181) The transformed Hamiltonian 7"^a-aa has no excitations with transfer energies 
|£fc,A-AA — £fc',A-AA| > A — AA. Let us first consider the transformation of 7io,A 

oo 1 

e x ^H ,xe- x ^ - H ,x = ^ e M ^ - { X^c+c* (Al) 

where we have introduced a new superoperator Xa,aa which is defined by the commutator 
of the generator X x< ax with operators A on which X Ai aa is applied, X A) aa-4 = [-Xa,aa> -Al- 
Note that the operators c\c k are the only part of Ho,x which do not commute with Xa,aa- 
We successively evaluate the commutators on the right-hand side of flAip . By use of 

XA,AAclc fc = ~^J2A~k ®T (^"^ *"* 4 C * + eKk '~ k)R% C U) (A2) 

k,i 

we first obtain 

Xa,aaW 0) a = ^ Akk ' ^ A ( £ k',x - e k ,x) e l{k - k '^ 5n t c\c k , (A3) 

kk'i 

Next by applying Xa,aa on Xa,aaWo,a where products of two fermions and of two local 
density operators Srij 5m occur. In order to keep only operators which are also present in 
7~£{x-ax) a factorization approximation has to be performed. With 

Xaaa 6m c\c k , = I]£ A~ kk e* k f x e^~ k ^ ((5n t 5n 3 )clc k , (A4) 

+ 5 k\k Sn i 8n j (cl,c k ,) - 5 kl - k (5n i 5n j )(clc k >) S ) 

k,j 

+5 kk 5n i 5n j {c\c k ) - S^SmSrij)^^)) 

one finds 

Xa,aaWo,a = ^ Y^ A kk>®kv X ( £ k>,\ - e k ,x) cos[i(k' - k)^ - Rft (A5) 

kk'ij 

x M6ni8nj)c~Ck' + 5 k , k 5n i 5n j (cl l c k ') - 8 k , tk -(5ni5nj)(c~c k >)) 
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or by using translational invariance 

x a,aa^o,a = ^ A 2 kk ,Q x k 'y X (e k ^x - e k ,x) (A6) 

kk' 

x (c p (k' - k)c\c k + (clc k )-^J2 5n * 6n i e^-W-^ - {c\c k )C p {k' - fc)) 

Note that in deriving (1A5|) and (1A6I) it was assumed that the number of k points which are 
integrated out in the step from A to (A — AA) is small compared to the total number of k 
points. Therefore, the product of two product 5 functions Q k k' and 6^, is only nonzero 
when the wavevectors are pairwise equal to each other (shell condition). In the next step, 
Xa,aa is again applied on (^ k c k . Thus by use of (IA2I) one finds 

Y 3 n, _ 1 r 2 2 C p (k - k')A 2 kk , , A AA j(k-k')R, A t 

x a,aa^o,a - I Jj \ Akk> ®kk> [ £ k',x ~ £k,x) e dUi C k C k ' 

kk'i 

(A7) 

which has again the structure of Xa,aa7^o,a except of the additional prefactor in the bracket 
[...]. thus one may sum up all commutator terms from (jAll) in order to obtain a compact 
expression for (lAlj) . In the same way one may proceed with the transformation of TC^x- 



APPENDIX B: GREEN FUNCTION Q Q {u) OF THE LATTICE GAS MODEL 

To evaluate the Green function Qo(uj) 

POO 

QoH = i / {With*)*, (Bi) 

Jo s 
we best use the Mori-Zwanzig projection formalism. First, let us define an anticommutator 
scalar product for operators A, B by (*4.|£>) = \[A^ , ■ 111 addition, we introduce the 

Liouville operator Lj belonging to the lattice gas Hamiltonian Tif. It acts on any operator 
A of the unitary space as L/A = [Hf,A\. Thus (IB II) can be rewritten as 

POO _ 1 

=i dt e* 1 *** {f}\e-^ * //) = (//| ? fl) (B2) 

JO Lf — UJ — 17] 

The well-known operator identity of the projection formalism 20 can be used to transform 
Qo(uj) into 

U — UJ — 2->\0j) 
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where the frequency term Q and the selfenergy E(u;) are given by 

n = (//£///) = (L///|g- r — QL///) (B4) 

Note that in (pj) and (El) we have already used that = = 1. The 



quantity Q is like Lj a superoperator which acts on usual operators of the unitary space. Q 
is a projection operator and projects on the subspace of the Liouville space perpendicular 
to |//), i.e., it is defined by Q\A) = \A) - \f})(f!\A). To evaluate Q and S (u) we need 

L/l//) = + 2 £ &|&i*/f) and QL/|//) = 2 £ ^QK//) (B5) 

where in the second relation the first term of L/|//) drops due to the presence of Q. As is 
easily seen, we obtain Q — if and 

E(w) =4££& gy (^//IQ-^— Q|^//) (B6) 

Note that the local density operators 5rii and 5rij in the bra- and ket vector of (1B6|) commute 
with Tif. Therefore, we can approximate £(a>) by an expression proportional to Qq 

= * 2 = 4 E£^ & ( Sn i Sn i)H f ( B7 ) 

which is identical to (jUJ). Note that in flB70 is independent of the lattice site i. £(u;) 
together with (1B3|) forms a quadratic equation for Qq{oj) which can be solved. It is easily 
seen that the imaginary part is given by 



for the frequency range \uj — if\ < 2k and Qq{ui) = elsewhere. Thus, 5sQo(u)) is nonzero 
only in a small frequency region around uo = if of width ±2k where k is determined by 
static density- density correlations of the lattice gas model. The maximum value of ^sQq(uj) 
is located at u = if and is given by 1/k. 

From ^sQ (u) the functions Bq(u) and Bq(u) are determined by help of the fluctuation- 
dissipation theorem 

B+(u) = (1 - f(u))-Qg (u) B (lu) = f(u)-Qg (u) (B9) 

7T 7T 
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where f(u)) is again the Fermi function. From (1B9|) and (j4T)l) . (j4"I]) the /-electron Green 
functions of the original Falicov-Kimball model can be found. 
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